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In this paper the general integral equations governing the mode spectra 
of optical masers are investigated from a point -of -view based upon certain 
theoretical results for Holder continuous kernels. Using an estimation origi- 
nally performed by Fredholm, it is proved that the homogeneous integral equa- 
tion 



*{x) = A f K(x,y) <p(y) dy 



has at least one eigenvalue for Holder continuous kernels K with exponent 
a > £ and with nonvanishing trace. A 11 the integral equations which have 
been treated in laser theory so far can be "factored" into one-dimensional 
equations with continuously differ enliable kernels, to which this result ap- 
plies directly. 

A Ithough in practice the vanishing of the trace is the exception rather than 
the rule, the later sections of this paper are devoted to demonstrations of the 
nonvanishing character of the trace of several of the common "laser kernels" 
associated with practical reflector configurations. These results provide in 
almost all cases the first rigorous proofs of the existence of eigenvalues and 
eigenfunctions for the integral equations of the optical maser. 

I. INTRODUCTION 

Homogeneous linear Fredholm integral equations with nonsingular 
kernels of normal type (which includes Hermitian kernels as a special 
case), i.e., kernels for which 

[ K(x,z) R(y,z) dz = [ K(z,x) K(z,y) dz, 

have been extensively studied. Within the framework of complex- valued 
£ 2 functions, questions of existence, uniqueness, and representation of 
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solutions can be largely answered for such equations. Recently, however, 
a number of integral equations involving kernels which are neither 
Hermitian nor normal have arisen in laser theory. These kernels, of which 

K{x,y) = e mx -" )2 , (1) 

with k a given complex constant, may be considered representative, do 
have the seemingly beneficial property of being complex-symmetric, viz. 

K(x,y) = K(y,x). 

Unfortunately, due to the lack of a sufficient analytic theory for such 
kernels, this "advantage" has yet to be adequately exploited. Thus, 
although it is well-known that every Hermitian kernel distinct from the 
zero transformation has at least one eigenvalue, the existence of eigen- 
values for general complex-symmetric or other non-normal kernels 
still remains an open mathematical question. 

These remarks are not meant to imply that there has been a paucity 
of theoretical investigation to complement the widespread experimental 
work with masers and lasers of various geometries. Quite the contrary! 
Boyd and Gordon, 1 for instance, have shown that for the confocal geome- 
try the resultant integral equation is equivalent to that considered earlier 
by Slepian and Pollak and has prolate spheroidal wave functions as 
eigenfunctions. Somewhat later, Boyd and Kogelnik 3 generalized this 
work to resonators with unequal reflector apertures and curvatures. 
Moreover, iterative computational methods have been applied by Fox 
and Li 4,6 and Li 6 to integral equations arising from a wide range of inter- 
ferometer geometries. Their techniques have produced plausible numeri- 
cal descriptions of the characteristic low-order modes and eigenvalues 
for the configurations considered. 

Even with the contributions represented by the above papers, how- 
ever, there still remains a dearth of knowledge, in a mathematical sense, 
about the eigenfunctions and eigenvalues (if any) of the homogeneous 
integral equations encountered in the general theory of the optical 
maser. The nature of some of the mathematical questions yet to be 
answered in this area was considered in an early 1963 paper by S. P. 
Morgan. 7 Since that time some progress has been made regarding the 
existence of eigenvalues for certain "laser kernels." Newman and 
Morgan, 8 by means of lengthy Taylor series techniques, have proved 
that kernels of the form 

K(x,y) = G(x) F(xy) H(y) 

with rather general G, F, H and with nonvanishing trace possess at 
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least one nontrivial eigenvalue. Other recent work 9 ' 10 has centered 
around the use of the natural Hilbert-Schmidt expansions of "planar" 
and "near-confocal" kernels in terms of their singular systems [Ref. 11, 
p. 142 ff.]. 

In this paper we want to assume a somewhat different approach based 
upon certain theoretical results for Holder continuous kernels. We first 
prove (in Section III) the following 

Theorem: Let the kernel K{x,y) be Holder continuous in either variable, 
with exponent a > \, for a ^ x,y ^ b. Then if the trace of K does not 
vanish, the homogeneous integral equation 



<p(x) = X / K(x,y) <p(y) dy 

Ja 



has at least one eigenvalue. 

The essential step in the proof is an estimation of the coefficients in the 
classical series representation for the Fredholm determinant of the 
kernel K(x,y), an estimation originally carried out by Fredholm him- 
self. 12 

We next observe that all the integral equations which have actually 
been treated in laser theory so far can be "factored" into one-dimensional 
equations with continuously differentiable kernels, to which the above 
theorem applies directly. Although we expect that in practice the vanish- 
ing of the trace is the exception rather than the rule, we devote the 
latter sections of this paper to demonstrations of the nonvanishing char- 
acter of the trace of several of the common "laser kernels" associated 
with practical reflector configurations. These examples are indicative of 
the ease with which the existence of eigenvalues and eigenfunctions 
can be rigorously established for many of the one-dimensional kernels 
arising in the theory of the optical maser. 

II. MATHEMATICAL PRELIMINARIES 

In general we shall consider complex-valued kernels K(x,y) defined 
on the bounded real domain a ^ x,y ^ b. Thus, where limits of integra- 
tion on integrals are not specified, the integrations are to be performed 
over the interval [a,b]. We shall also assume that K belongs to the class 
£ 2 , i.e., 



- U J a \K{x,y) \ 2 dxdy~] 



norm K = \\ K \\ = \J j \K(x,y) | 8 dx dy < oo , 

and that K{x,y) is a square-sum mable function of y for each value of 
x and conversely. 
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Our notation for composite kernels shall be 

KL = [ K(x,z) L(z,y) dz. 
Iterated kernels will be denoted by 



K" = Kit"'* 

= [ K(x,z) K"-\z,y) dz v^2 
with K 1 = K(x,y). In the same manner 

trace K = tr(K) = f K(x,x) dx 
and 

K = tr(K') = [ K"{x,x) dx 

J a 

= [ [ K(x,z) K"\z,x) dzdx. 

Reference should be made to Smithies 11 for further definitions and 
standard theorems on integral equations as needed. 

Certain notions regarding the characterization of entire or integral 
functions will also be of value in our work. In particular, recall that the 
order n of entire f(z) = zl "= a> n z n is given by 

.. 71 log n , > 

H = ImiSUPi -rj, rr. (2) 

n ^» log (1/1 a n |) 

Other definitions and results may be found in texts such as Boas. 13 

The property of Holder continuity is indicative of the smoothness of 
a given function. Kernels for which there exist positive constants A and 
a such that either 

I K(x,y) - K(z,y) I < A I x — z \ " for all x,y,z in [a,b] 

or 

I K(x,y) — K(x,z) I < A \y — z\ a for all x,y,z in [a,b] 

are termed Holder continuous in x or y respectively with exponent a. 
If a = 1 the functions are said to satisfy a Lipschitz condition, and thus 
Holder continuity is occasionally designated Lip a . It should be noted 
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that continuously differentiable functions automatically satisfy Lip- 
schitz conditions with a = 1. 



III. THE MAIN THEOREM 

Theorem: Let the kernel K(x,y) be Holder continuous in either variable, 
with exponent a > \, for a ^ x,y ^ b. Then if the trace of K does not 
vanish, the homogeneous integral equation 

<p(x) = A [ K(x,y) <p(y) dy (3) 

has at least one eigenvalue. 

Proof: The eigenvalues of (3) are the zeros, if any, of the Fredholm 
determinant D(\) associated with the kernel K(x,y). The classical 
series representation 12 for this entire function D(X) is 





D(\) = £ dX 




(4) 


where do = 1 and 








<- l H r ff~ 


• ( k( Si > S2 >- '•>*>) danism 


■ ds y 


(p ^ 1) (5) 


with 








K Ai ,s 2 , • • • , sA 


= det (K(si,Sj)) 







K(s! ,si) K(s 1 ,s 2 ) ■ ■ ■ K(s! ,s y ) 
K(s 2 ,sO K(s 2 ,s 2 ) • • • K(s 2 ,s v ) 



(6) 



K(S V ,Sl) K(S„,S<i) ••• K(Sy,Sy) 

We want to determine the order m of D(\) under the above hypotheses 
on the kernel K(x,y). Let us assume, therefore, that K is uniformly 
Holder continuous with respect to the second variable, that is 

| K(x,y) - K(x,z) \ <A\y - z\ a (7) 

with a > \. 

To estimate the coefficients d v ,* we first transform the determinant 
in (6) by subtracting the second column from the first, the third column 
from the second, etc., thus obtaining 

* This estimation was originally performed by Fredholm in 1903. u 
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det (K(Si ,Sj)) = [(si - s 2 )(s 2 - s 3 ) • • • (s„_i - s„)] a 



«n £12 

621 f22 



«■■! fi.2 



(i,,_i K(si , s v ) 

e 2 ,,_i K(Si,S v ) 
e,._i K(s v , S„) 



where 



m = 1,2, • 



K(s m ,s n ) - K(s m ,s n+i ) 

(S„ — S n +l) a 

, v; n = 1,2, ■ ■■ , v - \. 



Since | e m „ | < A by (7), and K(x,y) itself is bounded by continuity, 
a simple application of Hadamard's inequality yields 

| det (K( Si , Sj )) | < CV /2 1 ( Sl - s 2 )(s 2 - s 3 ) • • • (a^ - s„) | °, (8) 

where C is a constant. Inasmuch as the determinant of (6) is symmetric 
in the s,- , we may assume in deriving a further upper bound that 



b ^ «i ^ s 2 ^ 



> s„ > a. 



The right side of (8) is then maximized by spacing the Si uniformly 
between a and b. As a consequence we obtain 



det (K(sify)) ! < Cv 



iv_ v/2 



e^r 



(r> i) 



from which it follows that the estimate 

\d,\ < (const.)' v~' {a+i) 



(9) 



is valid for the coefficients of the power series of (4). 

The relation (9) implies that the order of the entire function Z)(X) 
satisfies 



v < 



which becomes less than 1 for a > f . Since 

.6 



d 



x ( K( Sl , Sl )d Sl = -tr(K) 



(10) 



* We have used Stirling's expansion for the factorial function. 
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does not vanish by hypothesis, D(\) must be a nonconstant entire func- 
tion of order less than 1 and hence must have at least one zero (see 
Ref. 13, p. 22 ff). It follows then that the integral equation has at least 
one eigenvalue.* Q.E.D. 

For entiie functions of finite order a general product expansion follows 
from the Hadamard factorization theorem. In view of the above results, 
therefore, D(\) may be written as the canonical product 

D(\) =fl(l-±) (11) 



where Xi , X 2 , • ■ • are the zeros of D(X) arranged in order of increasing 
modulus. This expansion converges absolutely and uniformly on com- 
pacta. 

If we differentiate the two representations of Z)(X) given by (4) and 
(11), set X = 0, and make use of (10) we obtain 

tv(K) = Zr- < 12 >t 

K=l \ V 

Thus kernels Holder continuous in either variable with exponent a > \ 
are one of the overlapping categories of non degenerate kernels for which 
the expansion (11) and hence the relation (12) is valid. Other classes 
include 

{%) Hermitian kernels with only a finite number of eigenvalues of 
one sign or the other (Mercer's Theorem), and 

(n) composite kernels of the form K = K\K* (Lalesco's result; see 
Ref. 15). 

IV. ONE-DIMENSIONAL EQUATIONS FOR THE OPTICAL MASER 

Careful analysis of an idealized diffraction model for the optical maser 
results in the following coupled integral equations for typical field quan- 
tities, such as the current densities: 5 



*i(zi#i) = Xi / Kn(xi,yi\xi,y2) $i(a,ifk) dx-rfyi, 
Js, 

$2( £2,2/2) = X 2 / K21 (: r 2 , \)i ; x x , y\ )<Pi( Xi,yi) dxidyx, 



(13: 



* If the order n of the Fredholm determinant D(\) satisfies < n < 1 then we 
can conclude that there exists a countably infinite set of eigenvalues of the equa- 
tion (3) since entire functions of nonintegral order have an infinite set of zeros. 

t This relation is not generally valid for arbitrary £ 2 kernels. In fact it may 
be inferred from results of Salem 14 that there are continuous symmetric kernelB 
with denumerably many eigenvalues for which 25— 1 (*A») does not even exist. 
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where 

K n (xi f yi;xt, y 2 ) = K 2l (x 2 , y% ; xi , yi) 

= exp ik [(xi - x 2 f + (y x - y 2 f (14) 

+ hi(xi , y x ) + lh{x 2 , y 2 )\. 

In these relations, & and S 2 are the mirror surfaces, hi(xi , t/i) and 
h 2 (x 2 , y 2 ) represent their respective departures from parallel planes, 
and k is a dimensionless parameter containing the wavelength as well 
as various geometrical dimensions such as the average mirror separation. 

The solutions $1,2 of (13) are eigenf unctions which describe the field 
distributions at the reflectors of the possible normal modes of laser 
oscillation; Xi and X 2 are the corresponding eigenvalues that specify 
the loss and phase shift which a propagating wave undergoes between 
the reflecting surfaces. Note in particular that these coupled equations 
( 13) are single-transit relations; that is, they give the field at each mirror 
in terms of the reflected field at the other. They can, of course, be com- 
bined into a single integral equation with a composite kernel. 

In the derivation of the above relations, the active maser material 
is assumed to be infinite in extent, homogeneous and isotropic. Diffrac- 
tion effects at the edge of the reflecting surfaces are neglected. Moreover, 
the separation between the reflectors, as well as the radii of curvature 
of these surfaces, is taken to be large compared to typical transverse 
dimensions. 

Although the integral equations (13) are two-dimensional, the pre- 
ponderance of analytic work on this problem has been under the addi- 
tional assumption that the laser kernels (14) can be adequately approxi- 
mated by a product of functions of single variables. If the field quantities 
are correspondingly decomposable, namely if 

*(**) = X(x) Y(y) 

or 

*(r,0) = R(r) 9(0), 

the general problem may be reduced to consideration of integral equa- 
tions in one independent variable. These equations take the form 

<Pi(xi) = Xi / Kn(xi,x 2 ) (p 2 (x 2 ) dx 2 , 

Vt(x») = Xi / K 21 (x 2 ,xi) <pi(xi) dxi, 

•'iii 
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with the single-transit kernels K given, for instance, by 

K u (xi , x 2 ) = K 2 i(x 2 , xi) = expik [(a* - x 2 ) 2 + piOi) + PiMl (16) 

Since these kernels are in general continuously differentiable functions of 
their arguments, we can use the main theorem of the preceding section 
directly to show the existence of eigenvalues whenever the trace is non- 
vanishing. 

V. ANALYSIS OF SPECIAL REFLECTOR CONFIGURATIONS 

In this section we examine the kernels associated with several well- 
known practical interferometer geometries and establish the general 
nonvanishing character of the traces. 

5.1 Plane Reflectors 

For rectangular plane reflecting surfaces which are mirror images of 
each other, the integral equation of interest is 

<p(x) = X f K(x,x') v(x') dx' 

with kernel 

K(x,x') = exp [ik(x — x') 2 ]. 

It is a trivial matter to verify that the trace of this differentiable func- 
tion is 2, from which we infer the existence of eigenvalues of the above 
integral equation. In fact, for this kernel one can use parity arguments 
to show there are at least two eigenvalues for k 7^ 0. 

5.2 Circular Plane Reflectors 

When the plane reflecting surfaces have circular cross section the 
integral equation kernel becomes 

K(r,r') = J n (krr') r' exp [ik(r 2 + r ,2 )/2] r,r' in [0,1] 

where J„ is a Bessel function of the 1st kind and nth order. The integer 
index n is indicative of an angular variation e'" which has been sup- 
pressed. The trace of the kernel K is proportional to 

T n = f J n {kT)e ikT dr, 

Jo 

and thus the nonvanishing of T„ will imply the existence of normal 
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modes of oscillation for the system. Using an integral representation of 
Jo , for example, we obtain 

T = f I f e ikTt (l - t 2 y h e ikT dtd T 

•>0 T J—l 

= 1 f 1 (!_**)-* C e ikT{l+t) drdt 

fJ-l Jq 

i r 1 rjka+t) ,-1 

= s/-, (1 -' 2 >-'[ L TT^] d( - 

It is easy to see that for Im (k) ^ 0, for instance, the real part of the 
integrand is negative almost everywhere. Hence the trace of the kernel 
does not vanish in this particular case, and we may draw our conclusion 
as to the existence of eigenvalues. 

5.3 Other Reflector Configurations 

For certain kernels, of course, there is little ,to be learned from applica- 
tion of our results on Holder continuous functions. Such is the situation 
regarding the kernel 

K(x,x ) = e 

associated with mirror image reflectors of square cross section, each 
having the curvature of a sphere centered at the center of the other 
reflector. As noted previously, this particular kernel gives rise to eigen- 
values and eigenfunctions which may be expressed in terms of prolate 
spheroidal wave functions. 2 

At the same time, however, it is advantageous that the eigenvalue 
existence question can be easily settled for more general reflector geom- 
etries which do not exhibit as beneficial analytic properties as the 
confocal configuration. In particular, the kernel 

K(x,x') = exp ik[(x — x') 2 + p(x) + p(x')], 

which pertains to mirror image square reflectors with shape function 
proportional to p(x), has eigenvalues if 

T = f K{x,x) dx = f e 2ikp(x) dx ^ 0. 

The vanishing of T for practical geometries would certainly seem to be 
the exception rather than the rule. 
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5.4 Composite Kernels 

The above examples show how our main theorem can be used to 
establish simply yet rigorously the existence of eigenvalues and eigen- 
functions for the one-dimensional laser kernels arising when the reflectors 
are mirror images of each other. For more generalized configurations in 
which the reflecting surfaces may be of unequal size and curvature, the 
applicable kernels are of a composite nature (see Refs. 3 and 5). In 
view of this one might choose to reason from Lalesco's results on com- 
posite kernels mentioned earlier rather than from our main theorem. 
This would be an acceptable method of attack. However, since (12), 
relating the trace to the sum of reciprocal eigenvalues, is valid in both 
situations, a verification of the nonvanishing character of the kernel 
traces is needed in either case. As a last illustrative example we shall 
provide this verification for the integral equations associated with 
asymmetric spherical reflectors of arbitrary curvature. 

Let «i = — h , a 2 = — b 2 , Pi(xi) = <W and p2(x->) = /3.r 2 2 . The one- 
dimensional integral equations (15) then become appropriate for analysis 
of an idealized interferometer having two rectangular mirrors of unequal 
size and unequal curvatures. As usual, these two coupled equations (15) 
can be combined into a single integral equation for either <pi or <p 2 . 
Moreover, this new integral equation may then be split apart into two 
subsidiary equations according to whether the eigenfunction modes are 
even or odd. The kernels resulting from this division are given by 



^(^)-»£;{ < B s l i sa}-pw^+^ i +-) 

+ 2.r' 2 (l +|8) - 2x'x]}dx' 
and have traces 

(K.,) = 2 J"* {"' ( C0S 2 ,f f ) exp (»feP(l + a) 

v ' ' Jo J-bo It sin 2kxx 



(17) 



tr 



(18) 



+ x' 2 (l + j8) - x'x]} dx'dx. 



It is easy to show that at least one of these traces is different from zero 
for real k and arbitrary curvatures a, /3. 
Note first that 



tr(K e ) - tr(K ) = 2 [ * f * exp {2ik[x 2 (l + a) 

.'0 J b 2 

+ x'\\ + 0)1} dx'dx 



88 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1965 



= 4 [/ ' exp {2ik(l + a)x 2 )dx 1 
• \f ' exp {2i*(l + P)x*} dx'\ . 



Now neither of the two bracketed terms on the right-hand side van- 
ishes, since the real parts of these Fresnel integrals are positive for real 
k, a, jS. Thus, the difference between the two traces, and hence at least 
one of the traces itself, is different from zero [one suspects, of course, that 
both of the kernels (17) have non vanishing traces]. Although this argu- 
ment gives no measure of the loss to be expected with any individual 
eigenmode, it does show that normal modes of oscillation exist for this 
arbitrary asymmetric spherical configuration. 
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